Composition of quantum operations 
and products of random matrices 

Wojciech Roga^, Marek Smaczynski^ and Karol Zyczkowski^'^ 

""^Instytut Fizyki im. Smoluchowskiego, Uniwersytet Jagiellonski, PL-30-059 Krakow, 
Poland 

^Centrum Fizyki Teoretycznej , Polska Akademia Nauk, PL-02-668 Warszawa, 
Poland 



dedicated to the memory of Ryszard Zygadlo 

Email: <wojciech.roga@uj .edu.pl>, <log_marco@poczta. onet .pl>, 
<karol@tatry . if . uj . edu . pl> 

May 6, 2011 

Abstract: Spectral properties of evolution operators corresponding to random 
maps and quantized chaotic systems strongly interacting with an environment 
can be described by the ensemble of non-hermitian random matrices from the 
real Ginibre ensemble. We analyze evolution operators 5* = ^'s ^'i repre- 
senting the composition of s random maps and demonstrate that their complex 
eigenvalues are asymptotically described by the law of Burda et al. obtained 
for a product of s independent random complex Ginibre matrices. Numerical 
data support the conjecture that the same results are applicable to characterize 
the distribution of eigenvalues of the s-th power of a random Ginibre matrix. 
Squared singular values of ^ are shown to be described by the Fuss-Catalan 
distribution of order s. Results obtained for products of random Ginibre ma- 
trices are also capable to describe the s-step evolution operator for a model 
deterministic dynamical system - a generalized quantum baker map subjected 
to strong interaction with an environment. 



1 Introduction 



Under the assumption of classical chaos the corresponding unitary quantum evo- 
lution, representing dynamics of an isolated quantum system, can be described 
[3 H] by random unitary matrices of the circular ensembles of random matrices 
[3]. If the quantum system is not isolated, but it is coupled to an environment, 
its time evolution is not unitary. In such a case one describes the quantum 
state by a density operator p, which is Hcrmitian, p = p\ positive, p > 0, and 
normalized, Trp = 1. The time evolution of such a system can be described 
by master equations [4], or by quantum operations [5], which correspond to a 
stroboscopic picture and discrete dynamics. 
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A quantum operation is described by a superoperator ^, which acts on 
the space of density operators. Let N denotes the size of a density matrix 
p. Then the superoperator is represented by a matrix \1/ of size A''^. Such 
a matrix is in general not unitary, but it obeys a quantum analogue of the 
Frobenius-Perron theorem [5], so its spectrum is confined to the unit disk. 
Spectral properties of superoperators representing interacting quantum systems 
were investigated in [71 [H |9] and also analyzed in an NMR experiment [10] . 
Similar properties exhibit also non-unitary evolution operators analyzed earlier 
in context of quantum dissipative dynamics [TTJ [T^] . 

Under the condition of classical chaos and strong decoherence the spectral 
properties of one-step evolution operators of deterministic systems do coincide 
with these of random operations [6] and can be described [13] by the ensemble 
of non-hermitian Ginibre matrices. All entries of such a random matrix are 
independent Gaussian variables. Since a superoperator describing one-step evo- 
lution operator can be represented as a real matrix [TJ], we are going to apply 
random matrices of the real Ginibre ensemble [THJ [11] [T^ . 

The main aim of this work is to study spectra of evolution operators de- 
scribing compositions of random quantum operations. Furthermore, we analyze 
s-step evolution operators representing quantum systems periodically interact- 
ing with the environment and compare statistical properties of complex spectra 
with predictions of the theory of non-hermitian random matrices. 

This work is organized as follows. In section 2 we recall necessary defini- 
tions of relevant ensembles of random matrices and briefly review recent results 
concerning statistical properties of their products. Properties of random maps 
and their compositions are analyzed in section 3. The model deterministic dy- 
namical system - a variant of the baker map interacting with an environment is 
studied in section 4. 

2 Non-hermitian random matrices 
and their products 

Consider a random square matrix G of size A^ of the complex Ginibre ensemble 
[18], generated according to the probability density 



This assumption implies that each entry Gmn of the random matrix is an inde- 
pendent complex Gaussian variable of a fixed variance = £,'^/N, where ^ is a 
free parameter which sets the scale. Eigenvectors of a random matrix G from 
such an ensemble are distributed according to the Haar measure on the uni- 
tary group, while complex eigenvalues Zi are described by the joint probability 
distribution 




(1) 
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From this results one can evaluate the density of eigenvalues in the complex 
plane. The density is rotationally symmetric and is a function of the moduls 
r = |z| of an eigenvalue [T51 [T7] . 

-<»)^^^ (3) 

where T{s^x) denotes the incomplete Gamma function, 

T{s,x) = / t"-^ cxp{-t)dt. 

In the asymptotic limit of large matrix size N the level density becomes constant 
inside the disk of radius R ~ ^, and decays exponentially outside the disk. This 
fact, known as the circular law of Girko [19], is conveniently formulated under 
the normalization = 1/N so that ^ = 1, for which the spectrum of a random 
Ginibre matrix of a large dimension covers uniformly the unit disk. 

Several recent applications including multiplicative diffusion processes |20| . 
macroeconomic time series lattice gauge field theories [22 and chiral en- 
sembles of random matrices [23l [24] increased interest in statistical properties 
of products of non-hermitian random matrices 1251 . Let Y denote a product 
of s independent square random matrices of size N form the complex Ginibre 
ensemble, Y — G1G2 ■ ■ ■ Gs- The density of the spectrum of Y is rotationally 
invariant in the complex plane [26], 

STT 

Here = denotes the product of scale parameters of each of the 

random matrices. For simplicity we shall assume that all s random matrices are 
characterized by the same variance, so that (, ^ ^l- The radial density of the 
eigenvalues reads 

= for r<e (5) 

s 

Based on the exact results for the level density for the case s = 1 of a single 
random matrix |27| [171 [H] is was suggested by Burda et al [29] to describe the 
finite size effects for the spectral density of a product of s random matrices by 
the dollowing ansatz, which involves the complementary error function erfc(a;) = 



Pjv(r) = P{r) ierfc(g(r - O^iv) . (6) 

Here q is an adjustable parameter, which does not depend on the dimension 
N, and the above form was reported 29, 30 to describe well the data obtained 
numerically by diagonalization of products of random Ginibre matrices. 
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Another way to describe a non-hermitian operator A is to study its singular 
values. Their squares are equal to the eigenvalues of the positive matrix AA'' . 
To set the scale it is convenient to renormalize such a matrix and define 

W=^ (7) 

such that TrT4^ = 1. 

Let {Xi}, i = 1,...,N denote the non- negative eigenvalues of W. The 
normalization implies that their sum is equal to unity, so if we use a rescaled 
variable, Xi = NXi its average value is equal to unity, (x) = 1. 

If the matrix A is random square Gaussian matrix from the Ginibre ensemble 
the level density P{x) describing the Wishart matrix W given in ([7|) is asymp- 
totically (for a large matrix dimension N) described by the Marchenko-Pastur 
distribution 131], 



1 /4 

FCi{x) ^ —\ - -1 for 0<a;<4. (8) 

If A is obtained as a product of s independent square random Ginibre matrices, 
A = G1G2 ■ ■ ■ Gs, the level density of the Wishart-like matrix W is given by 
the Fuss-Catalan distribution of order s [5H 133] • The name of this distribution 
is related to the fact that its moments are equal to the Fuss-Catalan numbers, 
often studied in combinatorics [34 . An explicit form of the FC distribution of 
order two, 

, , [^(27 + 3V81^12a:)^-6^] 



127r 



X3 



(27 + 3V81-12a 



valid for x G [0, 27/4] , was derived first by Penson and Solomon 35 in context of 
construction of generalized coherent states. More recently this formula was used 
in [361137] to describe singular values of random quantum states, the construction 
of which involves the product of two random matrices. Treating the sequence of 
Fuss-Catalan numbers as given one can solve the Hausdorff moment problem and 
find the corresponding probability distributions FCs{x). They can be written 
down explicitly |38j and represented as a combination of the hypergeometric 
functions [39] of the type sFg-i of the same argument. For instance, the Fuss- 
Catalan distribution of order three reads 

FT (t) — 1 F I—-L i X- i 3. 2J-t\ — 1 p /'I 1 5. 3 5. JLrA 

31-^^ V27i-£c3/4 3-^ 2 V 12' 4' 12' 2' 4' 256'^/ 47r3;i/2 3-^ 2 V 6 ' 2 ' 6 ' 4 ' 4 ' 256 "^Z 

1 oi^n (A 3 13-5 1- ILx) HQ) 

32y27r2;i/<i 3-^ 2 V 12 ' 4 ' 12 ' 4 ' 2 ' 256'^/ ' ^ 

The support of the Fuss-Catalan distribution FCs{x) of order s is formed 
by an interval [0, (s + ly^^ /s"] [32l |36]. Although this distribution was shown 
first to describe asymptotic distribution of squared singular values of a product 
of s independent Ginibre matrices, A = G1G2 ■ ■ ■ Gs, the same law describes 
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asymptotically the distribution of squared singular values of the s-th power of 
a square Ginibre matrix, A = [40] . 

Analyzing discrete time evolution of Hermitian density matrices in terms of 
quantum maps one copes with evolution operators which can be represented 
by a real matrix [S]. In particular, superoperators associated with random 
quantum maps [5] can be described |13| by real random matrices of the Ginibre 
ensemble. These matrices can be formally defined by the distribution ^ applied 
in the space of real matrices, so the argument of the exponent can be written 
as TiAA^ . To generate a random matrix pertaining to this ensemble one takes 
N'^ independent random Gaussian variables of the same variance and forms out 
of them a non-symmetric square matrix. 

Statistical properties of the real Ginibre ensemble are more complicated to 
analyze [15j [16] than in the complex case. For instance, for the real Ginibre 
ensemble the joint probability distribution of eigenvalues depends explicitly on 
their imaginary parts, so the level density is not rotationally invariant. In 
fact, there exist an accumulation of eigenvalues along the real axis, which is 
compensated by the repulsion of complex eigenvalues in vicinity of the real 
axis [15l[T7]. However, the number of real eigenvalues of a real random Ginibre 
matrix of size N scales as -\/iV, so that the non-uniform features of the spectrum 
can be neglected in the asymptotic limit N oo. Also the distribution of 
singular values of products of s real Ginibre matrices can be described [37] 
by Fuss-Catalan distribution of order s, originally applied [52] for products of 
complex Ginibre matrices. 

3 Random operations and their compositions 

A quantum operation is a linear map, p — ^ p' = ^(p), which sends the set of the 
density matrices into itself, so it preserves positivity and trace of the input state 
p. Any quantum operation acting on a d-dimensional state can be described 
by a unitary evolution U applied to an extended system followed by the partial 
trace over the environment £, 

^ip)^TTs(uip®W){iy\)U^). (11) 

Here Iv) G "Hm denotes the initially pure state of the environment, which is 
assumed to be M dimensional, so the unitary matrix U of size dM acts on the 
composite Hilbert space Hd Hm- 

If we reshape a density matrix p of size N into a vector of length (P, a 
quantum operation can be represented by a x cP matrix called superoperator. 
It is convenient to represent a density operator p of size N by its Bloch vector, 

d^--i 

p=Y.a,Y ■ (12) 

1=0 

Here 7* denotes the set of — 1 generators of the group SU((i), which satisfy re- 
lations Tr (7*7-') = (5*-' , while 7° = 1/^/N. As any density matrix is hermitian. 
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p = p\ the components of the Bloch vector are real, G K for i = 0, . . . , — 1. 
Thus the action of a quantum operation can be described as an affine trans- 
formation on the Bloch vector a representing the quantum state, a! = Cd+ k, 
where C is a non-hermitian distortion matrix of order (P — 1, while k is a trans- 
lation vector of length (F — 1. Using the Bloch vector representation one writes 
the superoperator ^ as a real matrix, 



1 

K C 



(13) 



The spectrum of the non-symmetric matrix C belongs to the complex plane. 
Since the quantum operation preserves the trace of the density matrix, Trvl'(p) = 
Tr{p) = 1, the spectrum of the superoperator belongs to the unit disk. 

Assuming that the matrix U € U{dM) in pT|) is taken randomly with respect 
to the Haar measure one obtains a random quantum operation [B]. In such a 
case the spectrum of the associated one-step evolution operator ^ was shown 
to consist of a single eigenvalue equal to unity, corresponding to the unique 
invariant state, and the remaining part localized in the disk of radius R — 1/ ^/M 
centered at the origin of the complex plane |13j . This characterization becomes 
exact for a large system size d. Since the matrix C is real, the spectrum of ^' 
is symmetric with respect to the real axis. As in the case of the real Ginibre 
ensemble [151 HB H] there exists a clustering of eigenvalues along the real axis, 
but this effect vanishes in the limit d — > oo. 

The number M, which determines the size of the disk of complex eigenvalues, 
equal to the dimension of the auxiliary subsystem, can be thus considered as a 
control parameter of the model describing the interaction of the principal system 
with the environment. Technically, M determines the rank of the Hermitian 
dynamical matrix |5|, which describes the quantum operation. 



3.1 Spectral density of a random superoperator 

We analyzed spectra of evolution operators corresponding to compositions of 
random maps. Let ^' = ^'s • • • ^2 • ^'i, where all s random maps are assumed 
to be independent. Figure [T] shows the spectra of such operators for s — 1,2 
and s — 3 for maps acting on a quantum system of size d = 20. To show the 
structure of the spectrum we magnified the scale accordingly, letting the leading 
eigenvalue z = 1 to remain outside the plot. According to the prediction of the 
Ginibre ensemble for s = 1 the distribution of the spectrum is close to be uniform 
in the disk of radius R — 1/VM (apart of the clustering of eigenvalues along 
the real axis), while its structure changes for larger s. 

Consider first the simplest case, s = 2, in which two random operations 
(fTT|) act successively. Assume that the first random operation ^I^i is due to an 
interaction with the environment £1 of dimension Mi , while the second operation 
\E'2 describes the interaction with the environment £2 of dimension M2. The 
resulting dynamics takes place in a tri-partite system described in the Hilbert 
space H = Hp Hi (g) 7^2- The first label p refers to the principal system 
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) s=2 C) s=3 

10-2 Re(A) xlQ-" 




-4 -2 2 4 xlO-2 -10 1 



Figure 1: Superimposed spectra of 25 superoperators C of dimension cP — 1 
associated with random maps acting on density operators of dimension d — 2Q 
and obtained by an interaction with an environment of size M = 20. The 
superoperators represent a) single random maps, s = 1; the composition of b) 
s — 2 and c) s = 3 random maps. The disk of radius Rs — (note 
the rescaUng of both axes) denotes the support predicted for the ensemble of 
random Ginibre matrices. 



of dimension d, while the other subsystems are labeled by the number of the 
environment £i and £2- 

After both operations the output state of the system, p" — '^2{p') — ^2[^i(p)] 
can be obtained by a three-subsystem unitary evolution, 

U = V2Vi = {Up2®lll){Upl®2l2). (14) 

followed by the partial trace over the collective environment £12 of dimension 
M = M1M2. Note that the resulting unitary U does not have a product struc- 
ture, as the symbols ® represent the tensor product with respect to two the 
different splitting of the Hilbert space, % — Hpi (82 ^-2 = ^1 T^i- 

It is natural to expect that in the case of two random operations, ^' — -^i, 
the resulting effect will be similar as the one caused by a single interaction with 
the combined environment £12 of dimension M = MiM2- This statement is 
equivalent to an assumption that the effect of the action of the resulting unitary 
U of the structure (fT4|) is statistically indistinguishable from the effect due to a 
global random unitary matrix Upi2 € U{dMiM2) followed by the partial trace 
over £12. Our numerical results confirm that this approximation works fine, as 
it implies that the subleading eigenvalues of the superoperator ^ live in the disk 
of radius R2 = 1/Vm = 1/^/M[M^. 

The same reasoning shows that in the case of a composition of s random 
operations, in which all dimensions of the environments are equal. Mi = M2 = 
■ ■ ■ = Ms = M, the essential spectrum of the resulting superoperator is asymp- 
totically confined in the disk of radius Rs = Rl — 1/Vm^. These predictions 
describe well the spectra of the evolution operators associated with the compo- 
sitions of two and three random maps and presented in Fig. [l]in panels b and 
c, respectively. 
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Figure 2: Radial density of complex eigenvalues of superoperators associated 
with a) a single random map, s = 1 and a composition of b) s = 2 and c) 
s = 3 random maps. Numerical data were obtained from a sample of 1000 
superoperators of dimension cf — 1 with d — 20. The dimensions of all auxiliary 
subspaces are equal, M = 20. Solid lines represent predictions ([S]) for radial 
density for products of s Ginibre matrices with the correction ([5]) due to finite 
size effects. Best fit gives the following values of the fitting parameter, a) g w 6 
for s = 1 , b) g w 7 for s = 2, c) g ~ 68 for s = 3. Dashed vertical line represent 
the radius Rg of the disk, which determines the support of the essential spectrum 
in the limit d — > oo. 

The radial distribution of complex eigenvalues collected of 1000 superim- 
posed spectra are presented in densities in Fig. [21 Numerical results can be 
described by the law ([5]) obtained in 26 for a product of s independent random 
complex Ginibre matrices. To take into account the finite size effects we used 
the ansatz ^ and fitted the parameter q. 

3.2 Singular values of random superoperator 

The distribution of squared singular values of a superoperator associated with 
a single random map or the composition of s of them was analyzed numerically. 
The data obtained presented in Fig. [3] show a fair agreement with the Fuss- 
Catalan distribution of order s, which describes properties of a product of s 
complex Ginibre matrices [32] . 

3.3 Average entropies 

To characterize the eigenvalue distribution P{\) of the positive matrix ^t'^J/t/XrvI'vl/t 

one studies the Shannon entropy of the spectrum, S = — A^lnAi. We 

put aside the leading eigenvalue Ai = 1 and renormalize the remaining cP — 1 
eigenvalues so that their sum is set to unity. The mean entropy computed 
numerically for a sample of 1000 superoperators 4* representing random maps 
with parameters d = M = 20 reads (S)^ ~ —0.505. This value agrees with 
the asymptotic prediction for the Wishart matrices, {S)i — implied by the 
Marchenko-Pastur distribution ([8]). A similar agreement is obtained in the case 
5 = 2, for which the average entropy implied by the Fuss-Catalan distribu- 
tion of order two © {S)2 = — | ~ —0.833 |36j, while the numerical data give 
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Figure 3: Density of normalized squared singular values, Xi = {cP — l)Ai, of 
superoperators corresponding to a) a single random map, s = 1; a composition 
of b) s = 2 and c) s = 3 random maps obtained for c? = 20 and M = 20. 
Numerical data collected from a sample of 1000 random maps are compared 
with the corresponding Fuss-Catalan distributions (|8I9[[T0)) represented by solid 
curves. 



the average entropy {S)<ii w —0.841. Numerical results for s = 3 provide the 
value {S)<^ K, —1.093, whereas the Fuss-Catalan distribution FC^{x) leads to 
(5)3 = -li « -1.083. 



Generalized quantum baker map 
and s— step evolution operators 



To investigate statistical properties of evolution operators associated with deter- 
ministic quantum systems interacting with an environment we shall concentrate 
on a model dynamical system, the classical analogues of which is known to be 
chaotic. Following the work of Balazs and Voros [41] we consider the unitary 
operator describing the one-step evolution of the quantum baker map. 



B = F] 



d/2 







Fd/2 



(15) 



Here Fd denotes the Fourier matrix of size d, namely [Fd\j-k = exp{\jk/2T:d)/Vd 
and it is assumed that the dimension d of the Hilbert space Hd is even. 

The standard quantum baker map B may be generalized, if in the defini- 
tion of the unitary operator (jlSp the Fourier matrix Fd is replaced by a two- 
parameter matrix p^'-''^'^^ 



\FS 



1 



Vd 



exp 



(j+(/)l)(fc-f (/)2) 



27rd 



(16) 



see Appendix D in [42]. The choice of both phases in [0,27r) does not influ- 
ence the classical limit, equal to the classical baker map. Thus combining Eq. 
(jl6p and (jlSp one obtains a two-parameter family of unitary quantum model 
dynamical systems, which we denote by B^^^^^- 

A certain variant of a non-unitary baker map introduced by Saraceno and 
Vallejos is capable to describe a dissipative quantum system [43]. Here we 
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are going to investigate yet another model of non-unitary quantum baker map 
introduced in [7l|44], which is deterministic, conserves the probabiUty, and is 
capable to describe projective measurements or a coupling with an external 
subsystem. Such a nonunitary dynamics can be represented as a quantum map 
and written in its Kraus form [S], 

M 

p' = 'i>{p) = Y,X,pX]. (17) 

For any trace preserving operation the set of M Kraus operators satisfies the 
identity resolution, X^jli -^l-^j ~ ^■ 

The parameter M is equal to the size of the environment coupled to the 
principal system of the baker map. Alternatively, M can be interpreted as the 
number of different outcomes of a measurement process. It is a free parame- 
ter of the model, which describes the degree of the decoherence in the system 
introduced by the nonunitary map (jl7p . 

The model of the classical baker map is chaotic and its dynamics can be 
characterized by the dynamical entropy of Kolmogorov and Sinai Hks equal to 
In 2 [3S] . To increase the degree of chaos one can simply take L iterations of the 
classical system, for which the dynamical entropy reads L In 2. This corresponds, 
in the quantum model, to taking the L-th power of the unitary evolution B 
given by (1151) . or its generalized version B^-^^^^ which involves (jl6|) . Increasing 
the parameter L one increases the degree of chaos in the classical model, and 
thus obtains unitary quantum operators which the properties of which are well 
described |46j by the Haar random matrices of the circular unitary ensemble. 

In our model system we take the i-step unitary dynamics of the baker map, 
B^ , followed by the non- unitary interaction with the environment of dimension 
M . In other words the non- unitary map $ given by (jl7p acts only every L steps 
of the unitary evolution. The structure of the complete evolution operator is 
presented schematically in Fig. |31 

Thus the stochastic quantum baker map describing the non-unitary evolution 
of the generalized quantum baker map [71 131] reads 

M 

It consists of M Kraus operators Pj, which act on the unitarily rotated state 
B'" p{B'^)^ . It is assumed that the ratio K = d/M is integer, so one can de- 
compose the Hilbert space T-Ld into the direct sum of M mutually orthogonal 
subspaces "Hq), j = 1, • ■ • , Af, of dimension K each. Then the Kraus operator 
Xj — PjB^ is a projection operator onto the fC-dimensional subspace 'H(j), 

so the sum X^jli -^l-^i ^ Sj=i -^o equal to identity, as required. Thus the 
parameter M in the non-unitary quantum baker map studied in this section has 
a similar meaning than M parameterizing the random maps: it describes the 
degree of the interaction of the principal system with the environment. Addi- 
tionally, for each choice of the system parameters (L, M) we may choose an pair 
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Figure 4: Sketch of the deterministic dynamical system - the generalized quan- 
tum baker map ^m.l ^ analyzed in this work: L steps of the unitary dynamics 
followed by an interaction with an M-dimensional environment £ described by 
the quantum operation $i. 
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of phases ((?!'i,^2) which enter (fTB)) and define the quantum model. To obtain 
a better statistics we shah superimpose spectra of superoperators obtained for 
fixed values of {L,M) and various phases {(pi,(j)2)- 

a) s=\ b) s=2 C ) s=3 



Rain> Re(A) Re(,\) 




-03 0.3 ^0.1 0.1 -4 -2 2 4 xlQ-^ 



Figure 5: Superimposed spectra of 60 superoperators of the generalized quantum 
baker map acting on density operators of dimension d — AO and characterized 
by parameters L — 20 and M — 10. The superoperators represent a) single 
baker map (fT8|) . s = 1; the s-step propagator {^m.lY for h) s — 2 and c) s — 3 
time steps. The disk of radius Rs = (note the rescaling of both axes) 

denotes the support predicted for the ensemble of random Ginibre matrices. 

Figure 5 presents exemplary spectra of superoperators of the generalized 
quantum baker map for s = 1,2,3 obtained for fixed values of L = 20 and 
M = 10 and 60 various pairs of phases (0i,02)- To display the structure of 
the bulk of the spectrum the scale increases with the power s, so the leading 
eigenvalue zi = 1 is located outside the figure. Apart of a few real eigenvalues, 
located for outside the circle of radius Rg = Af ^"Z^, the remaining eigenvalues 
are located close to the disk predicted for products of random matrices in the 
asymptotic limit d — >■ oo. The spectra are symmetric with respect to the real 
axis and exhibit the clustering of eigenvalues along the real axis combined with 
the repulsion of eigenvalues in the vicinity of the real axis. These effects, typical 
to the ensemble of real Ginibre matrices |15[|16) . vanish in the asymptotic limit, 
in which properties of products of complex and real random matrices tend to 
coincide. 

Figure 6 shows the radial density distribution P{r) for complex eigenvalues of 
superoperators of the generalized baker map for s = 1,2,3. The dynamical 
parameters of the model are fixed, d = 40, L — 20 and M — 10, while to 
accumulate a necessary statistics we superimposed data of 100 superoperators 
obtained for different values of the phases {(f>i, 02)- Note that already for s = 2 
the spectral properties of the superoperator $ can be described by the ensemble 
of random Ginibre matrices. 

To analyze properties of superoperators associated with s-step propagator 
of the generalized baker map, we analyzed also statistical properties of squared 
singular values of equal to eigenvalues of a positive operator Using 
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Figure 6: Radial density of complex eigenvalues of superoperators associated 
with a) one step evolution operator of the baker map, s — 1 and the s-step 
propagators for b) s = 2 and c) s = 3. Numerical data were obtained from a 
sample of 100 superoperators of the generalized quantum baker map acting on 
density operators of dimension d — AO and characterized by parameters L — 20 
and M = 10. Solid lines represent predictions ([S]) for radial density for products 
of s Ginibre matrices with the correction © due to finite size effects. Best fit 
gives the following values of the fitting parameter, a) q w 1.5 for s = 1, b) 
q w 1.5 for s = 2, c) g 3.5 for s = 3. Dashed vertical line represent the radius 
Rs of the disk, which determines the support of the essential spectrum in the 
limit d — > oo. 

the Kraus decomposition of a superoperator [5], $ = J^fLi-^i ® -^i ^^id 
that = X,X] (g, X,Xj. Since in our model each Kraus operator is 

a projector rotated by the same unitary matrix, Xi = PiU, the unitaries cancel 
out. The projections operators are mutually orthogonal, PiPj = SijPj so the 
above expression reduces to a single sum = X)j=i ^ Pj- ^'^Y operator 
Pi projects onto the subspace of dimension K — d/M so the operator = P, 
where P is a projector on a space of dimension N' = AI{d/M)'^ = d'^/M, so its 
spectrum consists of N' eigenvalues equal to unity and remaining d^{l — 1/Af) 
eigenvalues equal to zero - see Fig 7a. 

Consider now the case s = 2, in which we analyze the spectrum of $^($^)^. 
This operator can be written as $P$t — ($P)($P)t. Hence the singular val- 
ues of <i>^ are equal to the singular values of a truncated matrix $P. Thus the 
generic, non-zero eigenvalues of W2 = $^($^)^/Tr[$^($^)^] will be described by 
the Marchenko-Pastur distribution - see Fig. 7b. In a similar way, the singular 
values of are equal to the singular values of a truncated matrix c^^^^P so its 
squared singular values are described by the Fuss-Catalan distribution of order 
s' = s—1. Although the distributions FCs{x) are known to describe the asymp- 
totic distribution of squared singular values of a product of complex Ginibre 
matrices, they describe also statistical properties of (s — l)-step propagators of 
the generalized baker map. 
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Figure 7: Density of normalized positive squared singular values, Xi = {cP/M — 
l)Xi, of superoperators corresponding to a) one step evolution operator of the 
baker map, s = 1; the s-step propagators (<1'a/,l)'' for b) s = 2, c) s = 3, and 
d) s = 4 time steps (with the delta peak at a; = removed). Numerical data 
were obtained from a sample of 70 superoperators of the generalized quantum 
baker map acting on density operators of dimension d = 40 and characterized by 
parameters L = 20 and M = 10. Densities are compared with the corresponding 
Fuss-Catalan distributions of order s' = s — 1 (jH HI \TO\i represented by solid 
curves. 
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5 Concluding remarks 



We analyzed complex spectra of superoperators associated to compositions of 
s random maps. Statistical properties of the eigenvalues can be described by 
products of random matrices from the Ginibre ensemble. Due to the fact that 
a quantum map 5" preserves hermicity of a quantum state p, the superoperator 
^ can be described [T3] by an ensemble of real Ginibre matrices, for which 
a clustering of the eigenvalues along the real axis occurs [TS]. However, for 
large system sizes these finite-size effects can be neglected and the density of 
eigenvalues can be compared with predictions obtained for products of complex 
Ginibre matrices. In particular, the radial density P{r) of complex eigenvalues 
of the superoperator associated with the composition s random maps can be 
described by the algebraic law of Burda et al. ^Hl US]: while the distribution 
of the squared singular values of the superoperator (i.e. the eigenvalues of the 
positive matrix 'i>'$^ are described by the Fuss-Catalan distributions of order s 
[3211361 [311 Eg. 

Our numerical results support the conjecture that the distribution of eigen- 
values of a product of s independent random Ginibre matrices, G1G2 ■ ■ ■ Gg, 
obtained in [SO", 'SCT, describe also the spectrum of s-th power G^ of a given 
random Ginibre matrix G. This observation encouraged us to compare statisti- 
cal properties of s-step propagators of non-unitary quantum dynamical systems 
with the predictions of random matrices. Under the condition of strong classi- 
cal chaos and sufficiently large coupling with the environment the corresponding 
one-step evolution operators can be described by the ensemble of real Ginibre 
matrices [13]. 

Investigating a generalized version of a model dynamical system - the quan- 
tum baker map interacting with an environment - we demonstrate that statis- 
tical properties of complex eigenvalues of s-step evolution operators associated 
with such deterministic dynamical systems agree with predictions obtained for 
products of random Ginibre matrices. For the dynamical system investigated 
the operator ^^'^ is a projection operator with spectrum containing {0, 1}, so 
it cannot be described by random matrices. However, for a larger number s 
of the time steps, the squared singular values of can be described by the 
Fuss-Catalan distribution of order s' = s — 1 characteristic to the s'-th power 
of random matrices. 

Thus products of random non-hermitian matrices, used to describe matrix 
valued diffusion [5^ or random density operators [3S1 [37| , can also be applied 
to characterize statistical properties of multi-step evolution operators corre- 
sponding to generic quantum dynamical systems strongly interacting with an 
environment. 
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